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Abstract 
Motivated by the future precision test of the Higgs boson at an ete” Higgs factory, we calculate 
the production ete~ — ZHvy in the Standard Model with complete next-to-leading order elec- 
troweak corrections. We find that for \/s = 240 (350) GeV the cross section of this production is 


sizably reduced by the electroweak corrections, which is 1.03 (5.32) fb at leading order and 0.72 
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(4.79) fb at next-to-leading order. The transverse momentum distribution of the photon in the 


final states is also presented. 
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I. INTRODUCTION 


Recently, a Standard Model (SM)-like Higgs boson around 125 GeV was observed by 
ATLAS and CMS collaborations at the we 


, 2l. This discovery is a great step towards 
the understanding of electroweak symmetry breaking of the SM. So far, most measurements 
of the properties of this new boson are consistent with the SM prediction. The new physics 
that affects the Higgs couplings has been cornered to a decoupling region i Besides, 
since many extensions of the SM (like the supersymmetric models) contain a SM-like Higgs 
boson |5| whose properties can be quite similar to the SM Higgs boson, it is difficult for the 
LHC to verify whether or not this new boson is the SM one. In order to precisely study this 
newly discovered Higgs boson, an ete~ collider, the so-called Higgs factory, is needed. 

In such an ete” Higgs factory, the properties of Higgs boson can be measured with rather 
high precisions ld le). The dominant Higgs production is the Higgs-strahlung process e*e~ — 
ZH, where the ZH events can be inclusively detected by tagging a leptonic Z decay without 
the assumption of the Higgs decay mode. The individual Higgs decay branching ratios can 
then be directly measured as the fractions of the total ete” — ZH cross section by observing 
the specific states. For \/s ~ 240—250 GeV with an integrated luminosity of 500 fb~!, about 
O(10°) Higgs bosons can be produced per year, which allows to measure the Higgs couplings 
at a few percent |8|. So the electroweak radiative corrections should be taken into account 
in the theoretical calculations of the production rate. For the process ete” — ZH, the 
leading order calculation was performed in |9] and the one-loop electroweak corrections were 
calculated with the soft-photon approximation in [10H12] (a compact analytical formula for 
the electromagnetic corrections was given in |11] and a numerical calculation algorithm for 
the real photon emission was proposed in {13]). 

For an e+e” Higgs factory with \/s ~ 240—250 GeV another possibly important process is 
ete — HZ+. On one hand, it is an important part of the inclusive process ete” — ZH +X 
or can be distinguished for a hard photon; On the other hand, since the H Zy vertex occurs 
at one loop in the SM, the H Zy couplings is particularly sensitive to possible new physics 
contributions, such as the existence of new heavy particles propagating in the loop id ią. 
In this work we calculate this production in the SM with the complete next-to-leading order 
electroweak (NLO EW) corrections. In Sec. II we will give a description for the analytic 


calculations. The numerical results and discussions are given in Sec. III. Finally, we draw 


FIG. 1: The pentagon diagrams for the process ete” —> HZy. 


our conclusions in Sec. IV. 


II. A DESCRIPTION OF ANALYTICAL CALCULATIONS 


In the SM the process ete” — H Zy is induced by the electroweak interaction at lead- 
ing order (LO). Due to the small Yukawa couplings, we ignore the contributions from the 
Feynman diagrams involving the Yukawa couplings of light fermions. We denote the four- 


momenta of initial and final states in the process as 
e” (qu) +e (q2) > H (q3) + Z(q4) + (45) (1) 
The NLO EW corrections (Aggy) include two parts: 


e Virtual correction (Avir). 


We adopt the dimensional regularization to isolate the ultraviolet divergences (UV) in 


3 


the one-loop amplitudes. Then we remove the UV singularities by using the on-mass- 
shell renormalization scheme |16]. The pentagon Feynman diagrams in the calculation 
are presented in Fig[I] The reductions of N-point (N < 4) tensor integrals are imple- 
mented by using the Passarino-Veltman algorithm . But for the calculation of the 
5-point tensor functions, we adopt the Denner-Dittmaier method developed in Ref. 
to reduce the tensor integrals and use our fortran subroutines to perform numerical 
study, which has been validated in our previous works hd hd. We also numerically 
checked that our results are UV finite. 


Real photon radiation (Areal). 

Due to the exchange of virtual photon in the loops, the infrared (IR) divergences can 
appear in the virtual correction. According to the Kinoshita-Lee-Nauenberg (KLN) 
theorem , these IR divergences will be canceled by the real photon bremsstrahlung 
corrections in the soft photon limit. We denote the momenta of initial and final states 


for the real photon radiation process as 


e*(q1) +e (q2) > H(q3) + Z(qa) + ¥(G5) + 7(46). (2) 


We take the phase-space-slicing method bd to isolate the IR singularity in the 
above process. An arbitrarily small cut-off parameter 6, is introduced to split the 
phase space into soft region (Eg < 0s./s/2) and hard region (Eg > 65,/s/2). So the 


real photon emission correction can be decomposed into the soft and hard parts: 
Areal = Ad soft F AC hard. (3) 


In the soft photon approximation Bi we can calculate the soft part of the correction 


by using the following equation 


d8 z 2 
dAT soft = dovss i: d'd ( a č č ) (4) 
T? J g<ie/s/2 226 \N 96 3992 G6 


where Eg = ,/|q6|? + m3 and we give a small mass m, to the photon to eliminate the 
IR divergence (we checked that the dependence on this non-physical mass m, is exactly 
canceled when the real radiation correction and the virtual correction are combined). 
Since the hard part of the correction is insensitive to this fictitious photon mass, it 


can be directly evaluated by the numerical Monte Carlo method . We notice that 


4 


there are two photons in the real emission process and one of them should be tagged 
as the observed hard photon with pr > 10 GeV and |n| < 2. The phase space integral 


of these two identical photons in hard part of real emission can be expressed as: 


dq 2 
MIAE E 
~al ef se P6( Hy — Ey) 


+f Te | Meee za], (5) 


where the factor i is from the identical photons in the final states, and E, is the energy 
cut that corresponds to the above hard pr cut. In order to improve the numerical 
stability of Eq.(5), we adopt the method in the Ref. to carry out the integral 
Eq.(6). Since each of the two photons in the final states can be softer or harder than 
the other one with an equal probability, the Eq.(5) can be equivalent to: 


PB Ay 
I = 2x 6 
ji i J 2E; Aa 2B, 6) 


This means that we can technically assume the photon y(qs) to be the tagged hard 
photon and impose a transverse momentum “bd > 10 GeV and pseudo-rapidity cut 
26 


< 2 on 7(q5) in the numerical calculations 2 
n 


Finally, the total NLO EW correction of the process ete~ — H Zy is obtained by 
Adtot = Avir F AT soft F ACT hard- (7) 


We do the calculations by using the packages FeynArts-3.8 ls, FormCalc-8.2 and 
LoopTools-2.8 (we have the experience for using such packages Lidl bd lal ba) 32]). We 
analytically checked the gauge independence of our LO result by using the Ward identities. 
We also numerically checked our LO calculations in Feynman gauge(F.G.) with the package 
CompHEP 4.5.2 in unitary gauge(U.G.) [33]. In Tab.I, we present the comparison results and 
find that they are well consistent and are gauge-independent. In addition, by simultaneously 
interchanging the momenta and the polarization vectors of the two photons, we exploit Bose 
symmetry of the amplitudes of the real emission processes and found the values of the 
corresponding amplitudes do not change within numerical precision. We also numerically 
checked our result for the real radiation correction by using the Comphep program and 


found good agreement. 


250 2.172(2) 2.172(4) 
350 5.316(5) 5.316(9) 
500 3.562(3) 3.561(6) 
600 2.705(3) 2.705(4) 
800 1.708(2) 1.708(3) 
1000 1.184(1) 1.184(2) 


TABLE I: The comparison of our LO cross section of ete — H Zy in Feynman gauge with those 


calculated by CompHEP 4.5.2 in unitary gauge. 
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FIG. 2: The one-loop electroweak correction to the cross section of ete~ — HZy versus the 
soft cutoff logôs for My = 125.66 GeV and ys = 500 GeV: (a) showing respectively AChard, 


Avir + Adsoft and Aoyot; (b) showing Aoto¢ with the calculation errors. 
III. NUMERICAL RESULTS AND DISCUSSIONS 


In the numerical calculations we take the input parameters of the SM as 
m = 171.2 GeV, me = 0.519991 MeV, mz = 91.19 GeV, 
sin? Ow = 0.2228, a(m) = 127.918. (8) 


The Higgs mass is taken as my = 125.66 + 0.34 GeV lh, which is the combined result of 
the measurements of the ATLAS and CMS collaborations. 


We numerically check the stability of the results versus the soft photon cutoff parameter in 
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FIG. 3: The cross section of e'e~ + H Zy versus ys, showing respectively the LO result and the 
NLO EW corrections. The uncertainty caused by the 2c range of the Higgs mass (124.98 GeV 
< my < 126.44 GeV) is also shown (the shaded bands). 


Fig [2] where we assume ys = 500 GeV and m, = 1078 GeV. From the left panel of Fig Bit 
can be seen that the values of Avir, AChara and Ao sof, depend on the soft cutoff log ôs, while 
the total NLO EW correction Aio is independent of log ôs within reasonable calculation 
errors. Besides, we checked that the total correction is independent of m, for a fixed 0s. 
Therefore, in the following calculations we take the 6, = 2 x 107? and m, = 1078 GeV. 

In Fig] we plot the cross section of ete~ —> H Zy versus the center-of-mass energy vys, 
showing respectively the LO result and the NLO EW corrections. We can see that the 
production rate can reach a few fb in the threshold region \/s ~ 300 — 350 GeV (maximally 
it can reach 5.5 fb at LO and 4.8 fb at NLO), and the corresponding EW correction can reach 
—12%. For ys > 400 GeV, the cross section decreases rapidly due to the suppression of 1/s. 
At a 240 GeV Higgs factory, like the proposed LEP3 or China Higgs Factory (CHF), the 
cross section of ete~ + H Zy can reach 1.03 fb at LO and 0.72 fb at NLO (the corresponding 
EW correction is —30%), while at a 350 GeV Higgs factory, such as the ILC and TLEP, the 
cross section of ete~ — H Zy will reach 5.32 fb at LO and 4.79 fb at NLO (the corresponding 
EW correction is —9.8%). We can also find that the uncertainty of the cross section caused 


by the Higgs mass becomes small with the increase of \/s. 
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FIG. 4: The transverse momentum distribution of the photon at LO and NLO for the process 

ete” + HZy with ys = 240,350 GeV. The shaded bands correspond to the uncertainty caused 


by the 2c range of the Higgs mass (124.98 GeV < my < 126.44 GeV). 


Finally in Fig] we show the transverse momentum distribution of the photon in the 
process ete” — H Zy at LO and NLO for ys = 240,350 GeV. It can be seen that the NLO 
EW correction can greatly reduce the LO differential cross section at low pr region. The 
impact of the uncertainty of the Higgs mass on the pr distribution becomes weak as the 
collider energy increases. For \/s = 240 GeV most of the events are produced in the region 
of pp < 20 GeV due to the center-of-mass energy close to the production threshold; while 


for \/s = 350 GeV the pr value of the photon gets much harder. 


IV. CONCLUSION 


In this work we calculated the cross section of ete~ — ZHvy with complete next-to- 
leading order electroweak corrections in the SM. We found that for \/s = 240 (350) GeV the 
cross section of this production can reach 1.03 (5.32) fb at leading order and 0.72 (4.79) fb 
at next-to-leading order. In a future ete~ Higgs factory, this process can be measured as a 


precision test of the SM. 
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